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A procedure for the semi-blind separation 
and detection of digital co-channel signals, 
whereby after a simultaneous spatial-temporal 
equalization, iterative projections of the symbol 
vectors to be detected, in a vector subspace or 
a symbol alphabet O of finite size are initialized 
making use of a network user identifier IDi. 
These projections are then used to adapt the 
symbol vectors Sj to the vector subspace and 
to the symbol alphabet Q of finite size, making 
use of the method of least squares algorithm but 
performing only matrix multiplications without 
eigenvalue or singular decompositions. 
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1 

ITERATIVE PROJECTION WITH INITIALISATION USING A NETWORK USER IDENTIFIER 
Technical field 

5 

This invention relates to a procedure and a device for the semi-blind separation and 
detection of digital co-channel signals. Cellular mobile communication networks are 
usually limited by co-channel interference, i.e. the spatial reuse of the very same 
radio link becomes limited through co-channel interference. A radio link is defined 

10 by its frequency and/or its time slot (in the time-sharing multiplex method) and/or its 
code (in the code multiplex method). To allow more than one subscriber on the very 
same radio traffic channel, procedures have been proposed which are based on spa- 
tial separability and separation of the subscriber signals (SDMA, Space Division 
Multiple Access). For this purpose co-channel signals are purposefully generated 

15 and detected by the base station. 

Prior Art 

These procedures employ group antennas with connected signal processing, which 
20 by forming the antenna radiation pattern reduces the co-channel interference for in- 
dividual subscribers. Such a procedure is for instance described in US patent speci- 
fication 5,515,378 (Roy and Ottersten). Three fundamentally different methods can 
be distinguished. Those that are based on knowledge of the spatial structure of the 
antenna group (the so-called spatial-reference-procedures) as described in R. Roy 
25 and R. Kailath in "ESPRIT- Estimation of Signal Parameters via Rotational Invari- 
ance Techniques", IEEE Trans. Acoust., Speech, Signal Processing, Bd. 37, July 
1989, pp. 984-995; procedures that are based on the knowledge of the known signal 
sequence (so-called temporal reference procedure), explained by S. Ratnavel, A. 
Paulraj and A.G. Constantinides in "MMSE Space-Time Equalization for GSM 
30 Cellular Systems", Proc. Institute of electrical and Electronics engineers, IEEE, 
Vehicular Technology Conference 1996, VTC'96, Atlanta, Georgia, pp. 331-335; 
and so-called "blind" procedures which use known signal properties for separation 
and detection. 

35 In the case of spatial procedures, the angle of incidence of the separate signals by 
different subscribers is estimated with the help of extensive procedures, and the in- 
formation obtained is used to optimize the combination of the antenna signals. The- 
se spatially filtered signals are then detected in a separate receiver. Currently known 
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methods have a lower limit for the angular separability and fail below this limit. 
This is a difficult problem because ideally the lower limit should be zero degrees. 

A serious problem with the temporal reference procedure is the requirement to 
5 exactly synchronize the known training sequence with the incoming signals from a 
subscriber. Due to the multiple paths and the unknown distance between the 
transmitter and the receiver, the precise sampling instant is not known beforehand, 
and it also differs for the individual test signals. Only after synchronization can pro- 
cedures be applied to adapt the weighting parameters of the individual antenna ele- 
10 ments, which finally determines the effective radiation pattern. The alternative, the 
simultaneous synchronization and adapting of the weighting antenna parameters is 
computationally considerably extensive as shown in J. Fuhl, DJ. Cichon, E. Bonek, 
Optimum Antenna Topologies and Adaptation Strategies for SDMA, IEEE Global 
Communications Conference, Nov. 18-22, 1996, London, UK, pp. 575-580. 

15 

As far as the "blind" procedures are concerned, these use very common signal pro- 
perties such as transmitting with a fixed sampling rate (FSR), selecting the symbols 
to be transmitted from a finite alphabet size (finite alphabet, FA) and transmitting 
the signals with a constant envelope (constant modulus, CM), as explained in 

20 J.Laurila, E.Bonek, SDMA Using Blind Adaptation, ACTS Mobile Communication 
Summit, Aalborg, Denmark, October 7-10, 1997, pp. 314-319, or A-l van der Veen, 
A. Paulraj, Singular Value Analysis of Time-Space Equalization in the GSM Mobile 
System, Proc. IEEE ICASSP'96, May 1996, Atlanta, GA, pp. 1073-1076. These are 
computationally very extensive and have therefore not been used in the civilian 

25 mobile radio so far. The ILSP (Iterative Least Squares with Projection) procedure 
described in the publication by S.Talwar, M. Viberg, A. Paulraj called "Blind Sepa- 
ration of Synchronous Co-Channel Digital Signals Using an Antenna Array- Part 1: 
Algorithms". IEEE Trans. Signal Proc, Vol 44, pp. 1184-1197, May 1996, is based 
on the assumption of a finite number of transmission symbols, and is in principle 

30 suitable for the blind separation and detection of co-channel signals. It is, however, 
highly dependent on the accuracy of the initial values of the iterative signal estima- 
tion. The problem gets still worse when non-linear modulation methods are used for 
the transmission. This is unfortunately the case with modern mobile radio systems, 
for example in the GSM system (Global Systems for Mobile Communications) and 

35 in the DECT system (Digital Enhanced Cordless Telecommunications), which 
employ the non-linear Gaussian filtered Minimum Shift Keying modulation method. 
In such cases the above mentioned procedure must be complemented with an algo- 
rithm that takes advantage of the constant envelope of the signal, such as for instan- 
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ce ACM A (Analytical Constant Modulus Algorithm) proposed by A J. van der 
Veen, A Paulraj in "An Analytical Constant Modulus Algorithm", IEEE Trans, on 
Signal Proc, vol 44, pp. 1136-1155, May 1996. This requires further extensive 
arithmetic steps added to the blind detection and signal separation procedure, which 
5 in itself is already very complex. 

US patent specifications 5,619,533 (Dent) and US 5,519,727 (Okanoue) describe an 
equalizer which employs a MLSE (Maximum Likelihood Sequence Estimator) that 
obtains its information from a sole antenna signal. A space-time equalization is not 

10 possible as there is only one antenna signal, but on account of the enormous amount 
of memory needed for the intermediate steps of the MLSE, it would anyhow 
be disadvantageous. The patent specification US 5,432,816 (Gozzo) refers to 
a correction that together with the method of the least squares error provides re- 
cursive estimations of the channel. This equalizer method also only obtains its in- 

15 formation from a sole antenna signal. Receivers that employ such equalizer are 
based on inaccurate channel estimation values, because they for instance can- 
not utilize any structural properties of the transmitted signals. 

Disclosure of the invention 

20 

The invention is based on the requirement to provide a semi-blind procedure for the 
separation and detection of digital co-channel signals, comprising the steps: 

simultaneous (joint) spatial-temporal equalization, 

25 

initialization of iterative projections using the subscriber identifier, and 

carrying out the above mentioned projections for the adaptation of a symbol 
vector both to the signal subspace and a finite alphabet symbol constellation, in both 
30 cases making use of the least squares error algorithm, but performing only matrix 
multiplications without eigenvalue or singular value decompositions. 

The last step is repeated k times until the last iteration of the detected symbol vector 
no longer differs from the penultimate one. 

35 

The temporal synchronization for the initialization of the second step is easily 
achieved by correlating the respective subscriber identification with the basic vec- 
tors that span the desired vector subspace. 
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Another aspect of the invention is performing the projections in a decoupled way, 
dividing them between the individual subscribers. This allows the use of only simple 
matrix multiplication in these projections and avoids the excessive computing need- 
5 ed in traditional procedures. The extensive calculation of the pseudoinverse of the 
input matrix, including the orthonormal basic vectors of the desired vector subspace, 
is only required one single time, namely before initialization. 

The temporal equalization carried out simultaneously with the spatial equalization, 
10 which is based on the processing of several antenna signals, leads to a co-channel 
interference suppression distinctly improved from that which the temporal equaliza- 
tion of one antenna signal alone can provide. 

Short description of the drawings 

15 

The procedure is now more closely explained by means of the figures, in which 

Figure 1 shows the fundamental steps of the procedure and 

20 Figure 2 describes how to carry out the last two parts of the procedure in Figure 1, 
however with the generalization that the received signal matrix A (200) does not 
necessarily have to consist of orthonormal base vectors. The fact that the received 
signal matrix A (200) can also be used as an input makes the procedural step (60) 
necessary. 

25 

Preferred embodiments of the invention 

As shown in Figure 1, a matrix X of noisy input signals of several common channel 
signals serves as an input data record to be further processed by the procedure. 
30 The matrix X contains N sampling values, P times oversampled, of d co-channel 
signals that are derived from the M single elements of a group antenna. 

This received signal matrix X will be simultaneously spatially and temporally equa- 
lized in procedural step (10), in a known manner, so as to create a new matrix Y. In 
35 the preferred embodiment of the invention, the lines of the Y matrix and the lines of 
the symbol matrix S, having Block-Toeplitz-Form, span the same vector subspace. 
For this procedural step (10) either the singular value decomposition or the compu- 
tationally efficient but sub-optimal subspace tracking method can be used. This ex- 
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tensive calculation however is only required one single time in the course of the 
whole procedure. 

In procedural step (20) the initialization of the projection vector tj is carried out with 
5 the help of the network subscriber identifier IDj which is included in the correspond- 
ing (co-channel) signal S[. In every TDMA mobile radio system some known sym- 
bols for the channel estimation (e.g. midamble in the GSM system) are sent which 
can be used for the network subscriber identification. 

10 In the final procedural step (40) alternate projections of the matrix Y are done both 
on the vector subspace as well on the known symbol alphabet of finite size. Here the 
projections are nothing more than multiplication of this matrix with the projection 
vector tj and rounding of the symbol vector to the known symbols of the symbol 
alphabet £2 respectively. Additionally simultaneous adaptation of the symbol vector 

15 to the vector subspace and to the symbol quantity of finite size takes place by means 
of the least squares error algorithm. This last procedural step (40) provides after K 
iterations the vector including the detected symbols of the subscriber i. The iteration 
is stopped as soon as the symbol vector of two consecutive iterations is no longer 
different. One can then go on to the next subscriber i+1 and run the procedure again 

20 starting from step (20). 

Referring to Figure 2, the Moore-Penrose pseudoinverse A + is calculated in proce- 
dural step (50). 

25 If the matrix A was not formed by orthonormal base vectors using spatial-temporal 
equalization but the received signal matrix X was employed directly, the procedural 
step (60) becomes necessary to distribute the jammer and interference terms equally 
in space. This procedure is called "pre-whitening". 

30 In the known procedural step (20), the i:th signal is selected using a network 
subscriber identifier IDj, and the projection vector ti is initialized with it. In 
the initialization step (20) a substep (22) can be inserted for the choice of the opti- 
mal synchronization instant, which can be determined by the position of the maxi- 
mum cross correlation between the base vectors derived from the first procedural 

35 step (10) according to Figure 1 and the known network subscriber identifier IDi. 

In procedural step (70) the symbol vector corresponding to the i:th signal is compu- 
ted as the product of the projection vector and the received sfgnal matrix ti * A. This 
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symbol vector is projected on the vector subspace and on the finite set of the known 
symbols, making use of the method of the least squares error. Mathematically this 
procedural step can be calculated as Si^projaftj^A]. 

With this new k:th value of the symbol vector Sj (k) , a new calculation of the projecti- 
on vector t according to Si (k) =proj ft [ti (k-1) A] is carried out in procedural step (80). The 
iterations are continued until convergence is achieved, i.e. until an estimate of the 
symbol vector is no longer different from the preceding iteration. One can then pass 
on to the estimation of a next symbol vector corresponding to another signal. 

It should be noted that a special advantage of the procedure is that the number of 
signals to be so estimated is basically unlimited. However, only one single signal is 
detected at a time, i.e. the signal detection is always separate from the detection of 
all other signals. If the number of the subscriber signals in the cell is known, the 
procedure can be stopped after the estimation of these signals. Contrary to spatial re- 
ference procedures, it is unnecessary to estimate all signals in order to decide whet- 
her it is a usable signal or noise. One advantage of the procedure is in fact that the 
number of co-channel signals can be determined in the first procedural step in an 
elegant way. 

A further special advantage of the procedure is that it is independent of the employ- 
ed modulation format. It is also particularly suitable for the non-linear modulation 
methods which are usual in mobile communications of today. Such a modulation 
method does not, unlike known means, lead to an increase of the computing time. In 
addition, with our method it is possible to detect non-synchronized co-channel sig- 
nals. This procedure too is extremely robust against instantaneous changes of the 
interference sources, which means that no synchronization is required within the 
mobile communication network. 

It should also be noted that the procedure basically supplies two different signal esti- 
mates: The projection on the vector subspace, making use of the means of the least 
squares method, supplies a soft estimate (Soft Signal Estimation), while the result of 
the projection on the signal alphabet supplies a "hard" estimate (Hard Signal Esti- 
mation). The results of the soft estimation can be used for forward error correction. 

In real band-limited systems, oversampling in the time domain alone supplies too 
little information for blind separation of signals. One advantage of the present pro- 
cedure is the combination of signals from several or all antennas. This results in a 
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higher robustness against co-channel interference and moreover allows the estima- 
tion of all co-channel signals in FIR-MIMO (Finite Impulse Response, Multiple In- 
put, Multiple Output) applications. Further, our procedure allows a considerably 
simpler hardware structure than traditional beamforming methods. Therefore, a furt- 
5 her advantage consists in the fact that the estimation of the channel impulse res- 
ponse and the symbol detection can be carried out in one step. The channel estimati- 
on and the symbol detection usually occur separately. 

Unlike the traditional detection, using MLSE (Maximum Likelihood Sequence Esti- 
10 mation) together with the Viterbi algorithm, the signal separation in the present in- 
vention is based on iterative estimation. The projection of the estimated row space 
of the data matrix on the known, finite set of modulation symbols is used, employ- 
ing the means of the least squares algorithm. This procedure is computationally 
considerably simpler. 

15 

Traditional receivers employ only a short known part of the transmitted signal se- 
quences (so-called training sequences) to estimate the channel impulse response. In 
contrast to this, the data and channel estimation in our semi-blind procedure is car- 
ried out making use of all transmitted symbols. In this way a higher estimation accu- 
20 racy is achieved. 

A special advantage of this new procedure is the relatively short computing time for 
the semi-blind estimation, separation and detection of co-channel signals. Consider- 
ing the first part of procedural step (70), the multiplication of the projection vector t; 

25 with the matrix A requires 2 * $ * N arithmetic operations, when the matrix A is of 
the size 5 * N, with N being the number of sampling instants and 8 being the dimen- 
sion of the projection vector t. In general, the value of 5 is approximately the same 
as the number of the signals to be detected. The same number of operations are re- 
quired in arithmetic step (80) for the new calculation of the projection vector 

30 t. Because of this, K*d* (4*5*N) arithmetic steps are necessary, where d is the 
number of the estimated signals and K the number of the needed iterations per sig- 
nal. Extensive computer simulations have shown that already with K= 2.. .3 one 
finds convergence if known user identifiers and known basic vectors are used for 
initialization. The number N of sampling values required for the detection lies in the 

35 range of 40-90, which is about half of the symbols of one single GSM time slot. 
With these parameters our procedure, compared with spatial reference algorithms 
having the same bit error rate, requires only 1/3 of the antenna elements. 
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To further reduce the computational time of the embodiment according to Figure 2 
while increasing the robustness of the estimation of weak signals, the first symbol 
vector S! can after its final detection be subtracted from the received signal matrix 
A. This procedure is called serial interference cancellation. Here the best results are 
obtained if the received signals are estimated in the order of decreasing power. The 
improvement of the estimation robustness through serial interference cancellation is 
the best when the projection algorithm is applied directly on the data matrix without 
any estimation of the vector subspace. 

As a result, the co-channel interference is reduced for the detection of the remaining 
co-channel signals and the de-correlation step (50) is made less complex and 
the estimation of the subsequent symbol vectors is made more robust. 

In general, the here described procedures are not only computationally less exten- 
sive than traditional procedures, but they are also, using a suitable initialization 
based on subscriber identification, considerably more robust against noise interfe- 
rence and signals from neighbour cells. 

A further advantageous embodiment of the invention does not use any preceding 
estimation of the vector subspace at all. Normally a procedure for the separation and 
detection of co-channel signals would be extremely unreliable. However, because 
the initialization (20) with a subscriber identification supplies good initial values for 
the subsequent iterative projections, calculations can in this case too be performed 
with good results, which extensive computer simulations have shown. 

The means for the implementation of the described procedures is a signal processor. 
Signal processors, even if having small capability, are already used today in digital 
mobile communications receivers. For the appropriate implementation of the proce- 
dure in a signal processor, it has to be noted that the received signal matrix X is 
highly symmetrical. One column follows with time displacement the preceding one, 
while another element is added. This kind of symmetrical structure, which enables 
temporal equalization, also enables the block-by-block calculation of matrices 
and/or their inverse values in the first procedural step (10). Here the advantage is the 
largest, because this procedural step constitutes by far the largest part of computa- 
tion time through the necessary matrix inversions for the spatial temporal equaliza- 
tion. 
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In the procedural substep of the projection (70) of the symbol vectors Sj on the sym- 
bol alphabet £2, a further decrease in computational time can be achieved, as the ge- 
nerally complex elements of the symbol alphabet £2 are made real through derotation 
with differentially coded phase modulation. As a result, an implementation based on 
a signal processor that can only handle real value operations becomes possible. 
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Claims 

5 1 . A procedure for the semi-blind separation and detection of digital co-channel 
signals, characterized in that after a simultaneous spatial-temporal equalization, ite- 
rative projections of the symbol vectors Sj to be detected in a vector subspace and/or 
a symbol alphabet Q of finite size are initialized making use of a network user iden- 
tifier IDj, after which said projections are used to adapt said symbol vectors. Si to the 
10 vector subspace and to the symbol alphabet Q of finite size, in both cases making 
use of the least squares algorithm and performing only matrix multiplications with- 
out eigenvalue or singular value decompositions. 

2. A procedure according to claim 1, characterized in that the projections and 
15 their latest calculations are repeated until the detected symbol vector of the 

last iteration is no longer different from that of the penultimate one. 

3. A procedure according to claim 1 or claim 2, characterized in that the initiali- 
zation of every projection vector t\ is carried out for every co-channel signal separa- 

20 tely and decoupled from all other signals, employing a network user identity se- 
quence IDj for the respective co-channel signal. 

4. A procedure according to any one of the preceding claims, characterized in 
that the initialization of every projection vector tj is optimally synchronized when 

25 the position of the peak value of the cross correlation between the basic vectors, 
determined in the first procedural step (10) according to Figure 1, and the known 
network user IDj, is determined as the optimal synchronization instant. 

5. Procedures according to one or more of the preceding claims, characterized in 
30 that after successful detection, the detected symbol vector is subtracted from the 

received signal matrix X, Y or A in the order of decreasing signal power, estimated 
using signal matching filtering (correlation of the known training sequences with the 
received data), before the next symbol vector is detected. 
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